Extracting d-orbital occupancy from magnetic Compton scattering 

in bilayer manganites 
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We consider the shape of the magnetic Compton profile (MCP), Jmag{pz), in Lai.aSri.sM^Or 
for momentum transfer p z along the [110] direction and the associated reciprocal form factor B(r) 
defined by taking the one-dimensional Fourier transform of J mag {pz). B(r) is shown to contain a 
prominent dip at r fa 1 A, where the minimum value B min of B(r) can be related to the occupancies 
of the e 9 orbitals of d, x ?_ y i and d 3z 2_ r i symmetry in the system. We illustrate our procedure in 
detail by analyzing the measured MCP at 5K and the MCP computed within the framework of the 
local spin density approximation (LSDA) and comment on the differences between the measured 
and computed e g occupancies as a reflection of the limitations of the LSDA in treating electron 
correlation effects. 



I. INTRODUCTION 



The double-layered manganites La2-2xSri + 2 2: Mn207 
have drawn a great deal of attention as model systems 
that display a wide range of transport and magnetic prop- 
erties and undergo a variety of phase transitions as a 
function of temperature, doping and magnetic field. In 
Lai.2Sri.gMn207 (i.e., x = 0.40), investigated in this 
study, the CMR effect is a factor of ~ 200 at 129 K 
(just above the Curie temperature T c ~ 120 K) under 
a magnetic field of 7 T, and even at low fields the re- 
sistance changes by ~ 200%.i Neutron studied in this 
compound show that below T c the spins are aligned fer- 
romagnetically within the MnO planes, but are canted 
antiferromagnetically between adjacent Mn-0 planes. 

Compton scattering refers to inelastic x-ray scattering 
in the deeply inelastic regime and it is well known that 
this technique directly probes the ground state momen- 
tum density p(p) of the many body electronic system 
through a measurement of the so-called Compton pro- 
file, J(pz)i where p z is the momentum transferred in the 
scattering process^^ X-rays interact not only with the 
electron charge but also with the electron spin. There- 
fore, the majority (minority) spin Compton profile Jf 
(«7j_) can be extracted by measuring two Compton scat- 
tering cross sections: one with the direction of the mag- 
netization of the specimen parallel to the scattering vec- 
tor, and the other with the magnetization anti-parallel. 
Thus, in ferro- or ferrimagnetic materials one can zoom 
in on the magnetic properties via the magnetic Compton 



profile (MCP) J ma g(p) obtained by taking the difference 
between spin up (majority) and spin down (minority) 
Compton profiles 



JmagiPz) = J]{Pz) ~ J|(p«) • 



(1) 



Jrnag can be expressed in terms of a double integral of 
the spin density, pmagip) as 



J mag \Pz ) 



Pmag 

(p)dp x dp y 



(2) 



where p mag (p) = Pt(p) ~ Pl(p)- Tne possibility of using 
magnetic Compton scattering (MCS) to determine the 
momentum distribution of magnetic electrons was recog- 
nized quite early^ 5 ^ but since the scattering cross-section 
in the magnetic channel is typically several orders of mag- 
nitude smaller than for charge scattering, MCS experi- 
ments have become practical only in the last few years 
with the availability of high-energy, circularly polarized, 
X-rays at the synchrotron light sources^ 

In a recent MCS study of Lai^Sri.gM^Cv, we have fo- 
cused on the [110] MCP and shown that the shape of the 
[110] MCP contains a remarkable signature of the occu- 
pancy of the d x 2_ y 2 electrons^ By using a high magnetic 
field of 7T to maintain an electronically homogeneous 
ferromagnetic phase, we demonstrated that variations in 
the occupancy of the d x 2_ y 2 orbitals can be adduced from 
the [110] MCPs measured at different temperatures. The 
purpose of the present article is to expand on the discus- 
sion of Ref . concerning how the shape of the [110] MCP 
can be analyzed in terms of the reciprocal form factor 



2 



B(r), which is related to the MCP via a one-dimensional 
Fourier transform along the direction of the scattering 
vector. In this connection, detailed forms of the rele- 
vant B(r) functions for various t2 g and e g orbitals are 
given and their characteristic behavior with r is delin- 
eated. Since our purpose is to illustrate the fitting pro- 
cedure for interpreting the shape of the B(r), we only 
discuss the experimental B(r) at 5K and the correspond- 
ing theoretical curve obtained using the local spin density 
approximation (LSD A). We also comment on the discrep- 
ancies between these computed and measured B(r)'s as 
a reflection of the shortcomings of the LSDA in treating 
electron correlation effects. 

Concerning other MCP studies of the La-manganitc, 
Koizumi et alM^ have investigated the doping depen- 
dence of the [100] and [001] MCPs under a relatively 
low magnetic field of 2.5 T. They analyze their MCPs 
in terms of atomic and cluster type computations of the 
momentum density to gain insight into the occupation 
of magnetic orbitals in the system. All existing MCP 
data on the double-layered manganite have been taken 
at a momentum resolution of around 0.4 a.u. Recently, 
Mijnarends et alH have investigated the MCPs in La- 
manganite along the three high symmetry directions and 
have shown that the resolution of about 0.4 a.u. is not 
sufficient for investigating Fermi surface signatures in the 
MCPs. 

This article is organized as follows. Section [H] gives 
some technical details of our LSDA-based MCP compu- 
tations and briefly introduces the B(r) function which is 
the basis of our analysis. Section II I II presents and dis- 
cusses the details of how we analyze and fit the B{r) 
function. Section IIVI makes a few concluding remarks. 
As already indicated, in this article we only discuss the 5 
K measurements taken under a magnetic field of at 7T; 
for technical details of the experimental measurements, 
the reader is referred to Refs. andH. 



II. COMPUTATIONAL DETAILS 

In order to obtain momentum densities and Comp- 
ton profiles, the electronic structure of LaSr 2 Mn 2 07 was 
first obtained within an all-electron charge and spin self- 
consi stent KKR framework. The reader is referred to 
Refs. Il2ll3ll4ll5l for details of our band structure compu- 
tations. The formalism for computing m omentum densi- 
ties is discussed in Refs. Ilfill7ll8ll9l2'a The lattice data 
were taken from Seshadri et ali^ Exchange-correlation 
effects were incorporated within the local spin density 
approximation (LSDA) 22 . Doping effects were treated 
within the rigid band approximation where the band 
structure of La2-2a;Sri + 2 ;r Mn2 07 is assumed to be the 
same as that of LaSr2Mn207. More specifically, we have 
modeled the doped system by adjusting the Fermi en- 
ergy Ep in the majority spin band to accommodate the 
proper number of electrons for the x = 0.4 case^ 

A prominent feature of the band structure is the pres- 



ence of a large exchange splitting on Mn, which yields 
a nearly unoccupied minority spin Mn-3d bandiii*2iS4 
An analysis of the density of states shows that the ti g 
states lie at a binding energy of about 1 eV, while the e g 
states are placed essentially at the Fermi energy. There- 
fore our ab-initio calculation is consistent with molecular 
models which give a splitting of about 1 eV between the 
t2 g and e g levels^ We should keep in mind however that 
in a band structure scheme the t2 g and e g characters 
of Bloch wavefunctions vary continuously across various 
bands and that there is no strict gap between these levels 
for either the majority or the minority bands. 

If we assume that the 3.6 Mn electrons effectively re- 
tain their atomic character, then we would expect that 
these electrons will occupy the up-spin states following 
Hund's first rule. In particular, the first three electrons 
will go into the three lower lying up-spin ti g orbitals 
(d xy ,d xz ,d yz ), while the remaining 0.6 electrons will oc- 
cupy the e g orbitals (d x 2_ y 2, d 3z 2„ r 2). Moreover, since 
the t2 g and e g states are separated by a large energy gap 
of about 1 eV, which is much larger than fcT, the occu- 
pancy of t2 g states will change little with T, and the ef- 
fect of temperature will mainly be to redistribute the 0.6 
e g electrons between the d x 2_ y 2 and d 3z 2_ r 2 orbitals. A 
more quantitative handle on the occupancy of magnetic 
orbitals can be obtained in favorable cases by considering 
the reciprocal form factor, B(r), which is defined as the 
Fourier transform of the spin momentum density^: 



B(r) = J p mag (p)exp(-ip • r)dp . (3) 

We now express p m a g (p) as a sum over the momentum 
densities of the molecular orbitals ipf^ (p) of the mag- 
netic electrons, weighted by their occupancies nf. 

p™ ff (p) = E^i^ MO (p)i 2 - ( 4 ) 

i 

By transforming Eq. 0] into real space, it is straightfor- 
wardly shown that 

B(r) = / ^ MO (s)V\* M °(s + r)ds, (5) 

where the integrals on the right hand side give the au- 
tocorrelation of the magnetic orbitals, ipf' IO (r), i.e. the 
overlap between - l M (s) and the same orbital translated 
by a distance r. Alternatively, B(r) along a given direc- 
tion can be obtained directly from Eq. [3] by taking the 
ID-Fourier transform of the MCP along that direction. 
By modeling the right hand side of Eq. |5j and compar- 
ing the results to the experimentally determined B(r), 
insight into the occupation numbers can then be ob- 
tained. 
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FIG. 1: B(r) as a function of r along [110] obtained via a 
one-dimensional version of the Fourier transform (Eq. [3J of 
the experimental [110] MCP measured at 5K (thin line), and 
the corresponding computed MCP (thick line) after including 
the effect of experimental resolution broadening. 
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FIG. 2: Partial contributions B t 2 g , B x 2_ y 2 and B 3z 2_ r 2 of 
Eqs. I7l9l as a function of r. Note Bt2 g gives the total contri- 
bution from the three t^g orbitals of xy, yz, and zx symmetry. 
Computations are for Z e ff = 1. 



III. RESULTS AND DISCUSSION 

Figure n compares the B(r) along [110] obtained by 
using Eq. [3] from the measured MCP with the corre- 
sponding result from the computed MCP based on the 
LSDA. 26 Our focus is on understanding the nature of 
the pronounced dip around r w 1 A, which is seen in 
both theory and experiment. In order to interpret this 
minimum in terms of occupancy of magnetic orbitals, we 



model B(r) using Slater type orbitals (STOs).2L2£ The 
STO of xy symmetry for example is given by 

ip xy = N xye?q>(-£r) , (6) 

where £ is the exponent and N is the normalization fac- 
tor. Other STOs are of similar form with the factor of 
xy in Eq. replaced by one of a different symmetry. The 
overlap integral in Eq. [5] for various STOs is easily com- 
puted. The analytic form of the result is: 




where 

t = £r = (Z eff /3)r, (10) 

and Z e ff denotes the effective nuclear charge (all the 
quantities are implicitly in atomic units). Note that we 
have grouped the contributions of the three different t 2g 
orbitals (xy, yz, and zx) into a single quantity Bt 2g in 
Eq.Qfor convenience. By choosing an occupancy of 3 for 
the t 2g states, of / for the d x 2_ y 2 and (0.6 — /) for the 
d^ z 2_ r 2 state, the total B(r) is given by 

m , B t2g (r) + fB x 2_ y 2(r) + (0.6 - f)B 3z 2_ r2 (r) 
B{r) = 3~6 ' 

where we have imposed the condition B(0) = 1 since we 
compare MCPs normalized to one. 

Figures El and |21 provide insight into the modeling of 
B(r) based on Eqs. 6-11. Figure [3 shows the behavior 
of individual contributions B t2g , B x 2_ y 2, and £> 3z 2_ r 2, 
while Fig. [3] delineates the variations in the dip around 
r « 1 A as a function of the occupancy / of the d x 2_ y 2 
orbital. We have chosen Z e ff — 1 in all calculations in 
this article. It is clear from Eq.^|that Z e ff is a scaling 
parameter. Its precise value is not so important for our 
analysis since it merely serves to shift the position of the 
minimum in B(r), but does not affect the size of the dip 
in B(r) which is our main concern. The key result of our 
analysis is that the dip in B(r) in Fig. [21 which arises 
from the contributions in t 3 and i 4 , is produced mainly 
by the autocorrelation of the d x 2_ y 2 orbital. The d 3z 2_ r 2 
orbital is seen from Fig. [2Jto give a positive contribution 
near 1 A. With regard to the individual terms from the 
xy, yz and zx orbitals to B t2g (not shown in Fig. EJ, the 
d xy orbital gives a small positive contribution, while the 
d xz and d yz orbitals give negative contributions. 



4 



0.04- 



0.02- 



-0.02 ■ 



-0.04 ■ 



-0.06 ■ 




FIG. 3: Model B(r) function computed using Ea. llll is shown 
for various values of the parameter /, which gives the occu- 
pancy of the d x 2_ y 2 orbital, over the range / = 0.0 - 0.6 
(from top to bottom). Inset: Plot of minimum B m i n of B[r) 
as a function of / (see Eq. I12H . Diagram at bottom right: 
A schematic d x 2_ y 2 orbital centered at the origin with lobes 
pointed along [100] and [010] and the same orbital displaced 
along [110]. 



The minimum value B m i n of B(r) along [110] is seen 
from Fig. |3 to be correlated with the d x 2_ y 2 occupancy 
given by /. The inset to Fig. 0m fact shows that B m i n 
can be fitted linearly: 



B„ 



= af + 6, 



(12) 



where a = -0.130 and b = 0.004. The sketch of the 
d x 2_ y 2 orbital in Fig. is perhaps helpful in this regard: 
when the d x 2_ y 2 orbital is translated along [110], the pos- 
itive and negative lobes will overlap and yield a negative 
dip at a distance of the order of orbital dimensions. Fi- 
nally, we note that in writing Eq. |3| for B(r), we have 
not included the effect of finite experimental resolution 
on the MCP. The finite resolution Ap (FWHM) in mo- 
mentum space (assuming a gaussian resolution function) 
leads to a multiplicative attenuation of B{r) in r-space 
given by 



B a (r) =£(r)exp(-r 2 /S 2 ), 



(13) 



of Fig. d / = 0.35, which implies occupancies of 0.35 
and 0.25 for the d x 2_ y 2 and d 3z 2_ r 2 orbitals, respectively. 
The corresponding values for the 5K experimental B(r) 
of Fig.^are: 0.40 for d x 2_ y 2 and 0.20 for d^ z 2_ r 2w^ This 
preferential occupation of the d x 2_ y 2 orbital is somewhat 
surprising since the Mn-0 bond length for the apical O 
atoms is larger than that for the in-plane O atoms, and a 
simple molecular orbital scheme would suggest that the 
energy of the antibonding d 3z 2_ r 2 level will be lowered 
in relation to that of the d x 2_ y 2 level. Note, however, 
that electron correlations beyond the LSDA are needed to 
stabilize the canted antiferromagnetic order between the 
two Mn-0 planes observed in the absence of an external 
magnetic field. The occupancy of the d 3z 2_ r 2 orbital pro- 
motes ferromagnetism, 9 so that antiferromagnetic corre- 
lations will tend to reduce the d% z 2_ r 2 population. Ref- 
erence |30j argues that the on-site Coulomb correlation 
U of Mn d electrons can produce a significant splitting 
of the e g states, modify interlayer exchange interactions 
and magnetic ordering, and promote anisotropy in the 
electrical transport by reducing the conduction along the 
c-axis. Bearing all this in mind, we consider the discrep- 
ancy between the depth of the minimum between the 
LSDA and experiment in Fig.^to reflect the presence of 
electron correlation effects beyond the LSDA framework 
underlying our computations^*^ 

Turning to the behavior of B(r) beyond the dip around 
r ~ 1 A, Fig. ^ shows the presence of long range oscilla- 
tions in the experimental B(r) as well as the correspond- 
ing LSDA calculations (damping effect of the resolution 
notwithstanding). These oscillations at large r are ab- 
sent in our relatively simple STO-based model as seen 
from Fig. 3. We emphasize that sizable values of B(r) 
beyond atomic dimensions are a hallmark of electronic 
states extending over larger distances as a result of the 
mixing of Mn and O orbitals in the MnC>2 planes. These 
hybridization effects are crucially important for faithfully 
modeling the shape of the MCP in a global senseAiS Our 
relatively simple STO-based model is not meant to repro- 
duce such a global fit to the MCP, but to only highlight 
the atomic-like features hidden in the MCP via the use 
of the B(r) function at low r-values. As we have already 
pointed out, when the x-ray scattering vector lies along 
the [110] direction, the number of magnetic electrons of 
a specific symmetry, i.e., (i-electrons of x 2 — y 2 symmetry, 
yield a distinct atomic signature in the B(r), allowing us 
to extract the occupancies of the e g electrons. 



where 



S 



4y/hi(2) 
Ap 



(14) 



For Ap = 0.4 a.u. the exponential factor is 95% at r = 
1 A, therefore it has little effect on the value of B m i n . 

The minimum in the B(r) curves of Fig. ^ can now 
be used to derive the occupancies of the e g orbitals. Us- 
ing the fit of Eq. 1121 we thus obtain for the LSDA curve 



IV. SUMMARY AND CONCLUSIONS 

We discuss how the shape of the [110] MCP in 
Lai.2Sri.gMn207 can be analyzed in terms of the recipro- 
cal form factor B(r) given by the one-dimensional Fourier 
transform of the MCP. By interpreting B(r) as the au- 
tocorrelation function of the ground state wavefunction, 
one can obtain a handle on the occupancies of various 
magnetic orbitals in the system. In particular, the B(r) 
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along the [110] direction in La1.2Sr1.8Mn.2O7 contains a 
distinct dip at r w 1 A and the value B min of -B(r) at the 
minimum can be related directly to the occupancy / of 
the d x 2_ y 2 orbital as: B rn i n = af + b, where a = —0.130 
and b — 0.004. The specific forms of the B{r) functions 
associated with various tig and e g Slater-type orbitals are 
presented and used to analyze B(r) functions for the 5K 
MCP data and the computed MCP based on the LSDA 
as illustrative examples. In this way, for LSDA we ad- 
duce / = 0.35 or occupancies of 0.35 for d x 2_ y 2 and 0.25 
for d 3z 2_ r 2 orbitals; the corresponding occupancies de- 
rived from the experimental 5K MCP are 0.40 for d x 2_ y 2 
and 0.20 for the d 3z 2_ r 2 orbital. The larger experimen- 
tal value of the d x 2_ y 2 occupancy reflects limitations of 
the LSDA in accounting for electron correlation effects in 
the system. Our approach would allow a determination 
of the occupancies of magnetic orbitals through an analy- 



sis of the shapes of the MCPs more generally in complex 
materials. 
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